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ABSTRACT 

We present the discovery of seven ultracool brown dwarfs identified with the Wide- field Infrared Sur- 
vey Explorer (WISE). Near-infrared spectroscopy reveals deep absorption bands of H2O and CH4 that 
indicate all seven of the brown dwarfs have spectral types later than UGPS J072227. 51— 054031.2, the 
latest type T dwarf currently known. The spectrum of WISEP J182831.08+265037.8 is distinct in that 
the heights of the J- and H-h&nd peaks arc approximately equal in units of fx , so we identify it as the 
archetypal member of the Y spectral class. The spectra of at least two of the other brown dwarfs exhibit 
absorption on the blue wing of the ii-band peak that we tentatively ascribe to NH3. These spectral 
morphological changes provide a clear transition between the T dwarfs and the Y dwarfs. In order to 
produce a smooth near- infrared spectral sequence across the T/Y dwarf transition, we have reclassi- 
fied UGPS 0722-05 as the T9 spectral standard and tentatively assign WISEP J173835.52+273258.9 
as the Y0 spectral standard. In total, six of the seven new brown dwarfs are classified as Y dwarfs: 
four are classified as Y0, one is classified as Y0 (pec?), and WISEP J1828+2650 is classified as >Y0. 
We have also compared the spectra to the model atmospheres of Marley and Saumon and infer that 
the brown dwarfs have effective temperatures ranging from 300 K to 500 K, making them the coldest 
spectroscopically confirmed brown dwarfs known to date. 

Subject headings: infrared: stars — stars: low-mass, brown dwarfs — 

stars: individual (UGPS J072227.51-054031.2, WISEPC J014807.25-720258.8, 
WISEP J041022. 71+150248. 5. WISEPC J140518.40+553421.5, 
WISEP J154151. 65-225025. 2. WISEP J173835. 52+273258.9, 
WISEP J1828+2650, WISEPC J205628.90+145953.3) 



1. INTRODUCTION 

Brown dwarfs, objects with too little mass to sustain 
the high core temperatures necessary for stable ther- 
monuclear fusion of hydrogen, are the least massive, 
and possibly the most numerous products of star for- 
mation. ^AltjMiiglifirs_t 1 _rjre^ early 
1960's (IKumarl fl963l : IHavashi fc Nakanol fL963h . it was 
not until decades later that brown dwarfs were discov- 
ered in bulk by wide-area, red (700—1000 nm) and near- 
infrared (1—2.5 /Ltm) capable surveys such as the Two 
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Micron All Sky Survey (2MASS, iSkrutskie etld1|2006T) 
the Sloan Digital Sky Survey fSDSS. lYork et all 120001) 
and t he Deep Near-Infrare d Southern Sky Survey (DE- 
NIS, lEpchtein et all 119971) . The emergent spectra of 
brown dwarfs are so distinct from that of late-type M 
dwarfs that the creation of two new spectral classes, L 
and T (jKirkpatrick et all 119991: iBurgasser etall 120061) . 
was required in order to properly classify theni 11 !. The 
latest type T dwarfs currently known were discovered 
in the UKIRT Infra red Deep Sky Survey (UKIDSS, 
lLawrence et al.l I2007D and the Canada Fran ce Brown 
Dwarf Survey (CFBDS. [Delorme et al.ll2008bD and have 
estimate d effective temperatures (T e gs) of 500 to 700 
K (e.g.. IBurningham etal]|200l IDelorme et al~ll2008at 
iLeggett et al.l 120091: iLucas et al.ll2010HLiu et al.ll2011[) . 

Despite these spectacular successes, there exists a gap 
of nearly 400 K between the coolest spec troscopically 
confi rmed brown dwarfs at T e g ~ 500 K ( | Lucas et al.l 
120101 ) and Jupiter at T cff - 124 K (jHanel et all 119811 ). 
Although observations of star formation regions and 
you ng associations such as the Orion Nebula Clus- 
ter ([Weights et all |2009[ ). Chameleon I dLuhman et al.l 
120051), and TW Hydrae (jChauvin et al.l 120041 ) suggest 
that nature can form brown dwarfs that will eventu- 
ally cool to these temperatures once they have dispersed 
from their natal clusters, they are simply too faint to 

11 A comp endium of known L and T dwarfs can be found at 
http : //Dwarf Archives . org 
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have boon detected by the aforementioned surveys. Re- 
cently however, two brown dwarfs with estimated ef- 
fective temperatures of 300-400 K, WD 0806-661B 
(jLuhman et al.l[20TTI ) and CFBDSIR J145 829+101343B 
(CFBDSIR J1458+1013B lLiu et al.l l20"Tl. were discov- 
ered in targeted searches for companions to nearby stars. 
Although efforts to obtain spectra of these two com- 
mon proper motion sources have been hampered by ex- 
treme faintness in the case of WD 0806— 661 (J > 21.7; 
iRodriguez et al]l2011l ). and proximity to its primary star 
in the case of CFBDSIR J1458+1013B (<0'/l), their 
mere existence suggests that a population of nearby cold 
brown dwarfs awaits discovery. 

Foremost is the question of what these objects will 
look like spectroscopically and whe ther a new spec- 
tral c lass beyond T, d ubbed "Y" (jKirkpatrick et all 
Il999t iKirkpatrickl I2000T ). will be required in order to 
properly classify them. Chemical equilibrium calcu- 
lations and model atmospheres predict that as brown 
dwarfs cool below T e s ~ 600 K, their atmospheres pass 
through a series of chemical transitions which in turn im- 
pact t he appearance of t heir e mergent spectra (iLoddersI 
19991: iBurrows fe Sharp] 119991 : lLodders fc Feglevl l2002f 
Burrows et all 120031 ). At T eff ~ 600 K, the combina- 



tion and overtone bands of NH3 emerge in the near- 
infrarcct 1 !. At T c ff ~ 500 K, the prominent resonance 
absorption lines of Na I and K I in the red optical spec- 
tra of warmer brown dwarfs weaken as Na condenses 
out of the gas phase into Na2S and then K condenses 
into KC1. Finally, H2O and NH3 will also condense out 
at T c ff - 350 K and -200 K, respectively. Although 
each transition has been suggested as the trigger for 
the Y spectral class, focus has primarily been on de- 
tecting the NH3 bands because they are predicted to 
emerge at the hottest effective temperatures. Although 
NH3 absorption has been tentatively detected in the 
near-infrared spectrum of CFBD S J005910.90-011401 3 
(hereafter CFBDS J0059-0114, IDelorme et all[2008al) . 
this feature has not been confirmed in the spectrum 
of the coo ler object UGPS J 072227.51-054031.2 (UGPS 
0722-05: ITucas et al.l[2010l ). 

Independent of their spectral morphology, the study 
of these ultracool brown dwarfs will provide important 
insights into both ultracool atmospheric physics and 
the low-mass end of the stellar mass function. Be- 
cause brown dwarfs and exoplanets have similar atmo- 
spheric conditions, ultracool brown dwarfs are also ex- 
cellent exoplanet analogs that can be used as bench- 
marks for model atmospheres. The study of these ultra- 
cool brown dwarfs will therefore directly inform the in- 
terpretation and characterization of exoplanets detected 
with the next generation of hi gh-contrast imagers lik e the 
Gemini Planet Imager (GPI: iMacintosh et aLll2006l ) the 
Spectro-Polarimetric High-contrast Exoplanet Researc h 
instrument for the VLT (SPHERE; iBeuzit et ai]|2006D 
Proje ct 1640 at Palomar Observatory ( Hinklev et al.l 
1201 ID and the L- and M-b and Infrared Camera (LMIR- 
cam; iSkrutskie et al.ll2010( ) for the Large Bino cular Tele- 
scope Interferometer (LBTI). Simulations by iBurgasserl 

12 Although the fundamental band of NH3 has been detected 
in th e spectra of warmer T dwarfs at 10.5 fim (e.g., iRoellig et al.l 



2004; Cushing ct al. 2006), the intrinsically weaker near-infrared 
bands require a higher NH3 abundance, and thus lower effective 
temperature, to become a dominant opacity source. 



(|2004[) and lAllen et all ()2005[ ) have also shown that the 
space density of cold brown dwarfs is very sensitive to 
both the underlying mass function and the low-mass 
limit of star formation. Identifying and characteriz- 
ing a statistically robust sample of cold brown dwarfs 
will therefore provide two crit ical constraints on theories 
of low-mass star formation ([Whitworth fc Stamatelloll 
[2^0llWhit"worth et al.l 127)071) . 

One of the primary science go als of the Wide-field In- 
frared Survey Explorer (WISE, IWright et al] . a 
NASA mission that recently surveyed the entire sky at 
3.4 (Wl), 4.6 (W2), 12 (W3), and 22 (W4) /Ltm, is to 
identify such cold brown dwarfs. The Wl and W2 bands 
were designed specifically to sample the deep CH4 ab- 
sorption band centered at 3.3 /xm and the region rela- 
tively free of opacity centered at —4.7 /i m in the spec- 
tra o f cold brown dwarfs (see Figure 2 of iMainzer et al.1 
1 2 1 If ) . Since the peak of the Planck function at these 
low effective temperatures is in the mid-infrared, a large 
amount of flux emerges from the 4.7 /Ltm opacity hole, 
m aking the Wl— W2 colors extremely red (W l— W2 > 
2. IMainzer et al.l 120111 : iKirkpatrick et all 120111 ). Indeed 
such red colors are almost unique amongst astronomical 
sources making the identification of cool brown dwarfs 
with the Wl— W2 color alone relatively easy (see Figure 
12 of IWright et £0112010) . 

We have been conducting a search for cold brown 
dwarfs since the start of the WISE survey in mid Jan- 
uary 2010. This search has alread y born fruit with the 
discovery of six late-t ype T dwarfs (jMainzer et al.l [20 1 lb 
iBurgasser et al] 120111) two of which ha ve spectral types 
later than T8. IKirkpatrick et all (|2011[) present over 100 
new brown dwarfs, the vast majority of which have spec- 
tral types later than T6. In this paper, we focus on seven 
of the —100 brown dwarfs whose near-infrared spectra 
indicate that the they are the latest-type spectroscopi- 
cally confirmed brown dwarfs currently known. Indeed 
we identify six of these brown dwarfs as the first mem- 
bers of the Y spectral class. In §2, we briefly discuss 
our selection criteria before presenting the ground- and 
spaced-based imaging and spectroscopic followup obser- 
vations in §3. In §4 we present the properties of the 
first Y dwarfs, define the transition between the T se- 
quence and the Y dwarfs, and derive estimates of the 
atmospheric parameters of the new brown dwarfs using 
model atmospheres. 

2. CANDIDATE SELECTION 

The seven new brown dwarfs were identified as part of a 
larger ongoing search for cold brown dwarfs using WISE. 
A detailed descrip tion of this survey and ou r search crite- 
ria is presented bv IKirkpatrick et al.l (|2Qlll ). Briefly, can- 
didates were selected from the source working database 
derived from the first-pass processing operational coadds 
using color constraints derived from known T dwarfs and 
model atmospheres (in the case of brown dwarfs with T e ff 
< 500 K). Table □] lists the WISE designations and pho- 
tometry of the seven brown dwarfs, and Figure 1 shows 2' 
x 2' DSS I, 2MASS J and H, WISE Wl, W2, and W3, 
and W1W2W3 color composite images for each dwarf. 
Hereafter we abbreviate the numerical portions of the 
WISE designations as hhmmiddmm, where the suffix is 
the sexagesimal right ascension (hours and minutes) and 
declination (degrees and arcminutes) at J2000.0 equinox. 
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3. OBSERVATIONS 

The followup ground- and space-based observations of 
the seven WISE brown dwarfs are discussed in the fol- 
lowing sections. Although we present the near-infrared 
photometry of the brown dwarfs in this work for com- 
pleteness, we de f er the discussion of these data to 
IKirkpatrick et all ()2011h who present a more detailed 
discussion that places them in context with the larger 
population of brown dwarfs. In addition to the observa- 
tions of the seven WISE brown dwarfs, we also obtained 
a near-infrared spectrum of UGPS 0722—05 for compar- 
ison purposes. A log of the near-infrared photometric 
observations as well as the resulting photometry is given 
in Table [21 and a log of the spectroscopic observations is 
given in Table [3J 

3.1. Near-Infrared Imaging 
3.1.1. NE W FIRM /Blanco 

WISEP J1541-2250 was observed on the night of 2011 
Apr 17 (UT) with the NOAO Extremely Wide Field In- 
frared Imager (NEWFIRM) mounted on the Cerro Tololo 
Inter-American Observatory (CTIO) Victor M. Blanco 
4m Telescope. A description of the instrument, ob- 
serving strategy, and d ata reduction can be found in 
IKirkpatrick et all (|201 lh . The resultant J- and tf-band 
photometry is presented in Tabic [2] 

3.1.2. WIRC/Palomar 

Near-infrared images of WISEP J0410+1502, 
WISEPC J1405+5534, WISEP J1738+2732, and 
WISEPC J2056+1459 were obt ained using the W ide- 
field Infrared Camera (WIRC; IWilson et all l200l on 
the 200 inch Hale Telescope at Palomar Observatory. A 
description of the instrument, observing strategy, an d 
data reduction can be found in IKirkpatrick et al.l (|20 111 ). 
The magnitudes and / or limits for each brown dwarf are 
given in Table [21 

3.1.3. PANIC '/Magellan 

WISEPC J0148-7202 was observed on the night of 
2010 Aug 01 (UT) with the now decommissioned Pers- 
son's Auxiliary N asmyth Infrared Camera (PANIC; 
iMartini et al.ll2004T ) on the east Nasmyth platform at the 
Magellan 6.5 m Baade Telescope. A description of the 
instrument, observing strategy , and data reduction can 
be found in IKirkpatrick et al.l (|2011l) . The J- and H- 
band magnitudes of WISEPC J0148-7202 are given in 
Table H 

3.1.4. NIRC2/Keck II 

High resolution observations of WISEP J1828+2650 
and WISEPC J2056+1459 were obta ined with NIRC2 
behind the Keck II LGS-A O system (jWizinowich et all 
[2001 Ivan Dam etail l200l on the night of 2010 July 
1 (UT). A description of the instrument, observ- 
ing strategy, an d data reduction can be found in 
IKirkpatrick et al.l (|2011| ). The J- and if-band magni- 
tudes are given in Table [2] 

3.2. Near-Infrared Spectroscopy 
3.2.1. SpeX/IRTF 



A 0.9-2.5 /mi, low-resolution (R = A/AA « 150) 
spectrum of UGPS 0722—05 was obtained with SpcX 
(jRavner et al.l 12003ft on the 3 m NASA IRTF on 26 
Jan 2011 (UT). A description of the instrument, ob- 
serving strategy, and d ata reduction can be found in 
IKirkpatrick et all ([2011D . The spectrum, which is shown 
in Figured has a high S/N, reaching >50 at the peaks 
of the Y, J, and H bands. 

3.2.2. FIRE/Magellan 

Low-resolution (i?=250— 350), 1—2.4 /im spectra 
of WISEPC J0148-7202, WISEP J0410+1502, and 
WISEP J1541-2250 were obt ained with the Fold e d-port 
InfraRed Echellette (FIRE, ISimcoe et all l2(M 1201(1 
mounted at the auxiliary Nasmyth focus of the Magellan 
6.5 m Baade Telescope . A description of the instrument, 
observing strategy, and data reduction can be found in 
IKirkpatrick et all (|2011h . The spectra are shown in Fig- 
ure H 

3.2.3. NIRSPEC/Keck II 

WISEPC J2056+1459 was obser ved using the Near- 
Infrared Spectrometer fNIRSPEC; IMcLean et all 1 1 9981 
[2000] ) located on one of the Nasmyth platforms of the 10 
m Keck II telescope on Mauna Kea, Hawai'i. The 0'.'38- 
wide slit in the low resolution mode provides a resolving 
power of i?=2500. WISEPC J2056+1459 was observed 
with the N3 order sorting filter (1.143—1.375 /mi) on the 
night of 2010 Oct 21 (UT) and with the N5 order sort- 
ing filter (1.431-1.808 /mi) on the night of 2010 Nov 22 
(UT). 

A series of 300 s exposures was obtained at two dif- 
ferent positions along the 42" long slit. An AO V star 
was observed after each series of science exposures for 
telluric correction and flux calibration purposes. Cali- 
bration frames consisting of neon and argon arc lamps, 
dark frames, and flat- field lamps were also taken follow- 
ing the science exposures. The data were reduced in a 
standard fashion using the IDL- based REDSPE CFl re- 
duction package as described in IMcLean et al.1 ((2003). 
Since REDSPEC docs not produce uncertainty arrays, 
we generated them as follows. First we performed a sim- 
ple sum extraction using the rectified, pair-subtracted 
images generated by REDSPEC. We then scaled the 
spectra to a common flux level and computed the av- 
erage spectrum. The uncertainty at each wavelength is 
given by the standard error on the mean. The average 
spectrum is then corrected for telluric absorption and 
flux calibrated using the calibration spectrum generated 
by REDSPEC. Since the difference between the spectra 
produced by REDSPEC and our spectra was negligible, 
we used our spectrum for our analysis. Finally, the N3— 
and N5— band spectra were absolutely flux calibrated us- 
ing the WIRC photo metry (see Table [5]) as described in 
iCushing et all 00051 ) and merged to produce a 1.15 — 1.80 
/im spectrum. The final spectrum, which is shown in Fig- 
ure [21 has a peak S/N of 8 and 6 in the J and H bands, 
respectively. 

3.2.4. WFC3/Hubble Space Telescope 
13 See http://www2.keck.hawaii.edu/inst/nirspec/redspec 
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WISEPC J1405+5534, WISEP J1738+2732, and 
WISEP J1828+2650 were observed with the in- 
frared channel of the Wide Field Camera 3 (WFC3; 
iKimble et al.ll205§f ) on-board the Hubble Space Telescope 
(HST) as a part of a Cycle 18 program (GO-12330, 
PI=Kirkpatrick). The WFC3 uses a 1024 x 1024 
HgCdTe detector with a plate scale of 0'.'13 pixel -1 
which results in a held of view of 123x126 arcsecond. 
The G141 grism was used to perform slitless spec- 
troscopy of each brown dwarf covering the 1.07—1.70 
/im wavelength range at a resolving power of i?«130. 
For each brown dwarf, we first obtained four direct 
images through the F140W filter (A p =1392.3 nm) in 
the MULTIACCUM mode with the SPARS25 sampling 
sequence. Between each exposure the telescope was 
offset slightly. We then obtained four images with the 
G141 grism at the same positions as the direct images. 
The spectroscopic observations were also obtained in 
the MUTLIACCUM mode but using the SPARS50 
sequence. 

The raw images were first processed using the CAL- 
WFC3 pipeline (version 2.3) which not only subtracts 
the bias level and dark current but also flat fields the di- 
rect images (the grism images are flat fielded during the 
extraction process described below). Th e spectra were 
then extracted using the aXe software (|Kummel et al.l 
120091 ) . which is a suite of PyRAF/IRAF packages de- 
signed to extract spectra from the slitless modes of both 
WFC3 and the Advanced Camera for Surveys (ACS). 
aXe requires knowledge of both the position and bright- 
ness of the objects in the field of view. We therefore 
combined the four direc t images using MULTIDRIZZLE 
(jKoekemoer et al"1l2002| ) and the latest Instrument Dis- 
tortion Coefficient Table (IDCTAB). A catalog of the ob- 
j ects in the field was the n constructed using SExtractor 
(|Bertin fc Arnoutslll996D . For each object in the source 
catalog, two dimensional (2D) subimages centered on the 
first-order spectra of each object were then combined us- 
ing the task AXEDRIZZLE to produce a high S/N 2D 
spectral image. One-dimensional, flux-calibrated spec- 
tra and their associated uncertainties arc then extracted 
from the 2D drizzle subimages. 

Since the G141 grism mode is slitless, spectral con- 
tamination from nearby sources is not uncommon. The 
aXe software (using the Gaussian emission model) esti- 
mates the level of contamination for each object using 
the positions and magnitudes of all the objects in the 
field of view. The spectrum of one of the brown dwarfs, 
WISEP J 1828+2650, exhibits moderate contamination 
that increases in intensity towards shorter wavelengths 
(see Figure [3|) . The aXe software does not actually cor- 
rect for this contamination so we attempted to do so 
using the contamination image generated by aXe. Unfor- 
tunately, the contamination-corrected spectrum exhibits 
negative flux values which suggests that aXe is over esti- 
mating the contamination level. We will therefore use the 
contaminated spectrum and consider it an upper limit to 
the actual spectrum. This issue will be discussed in more 
detail in t j4T.ll 

4. ANALYSIS 
4.1. Spectral Characterization 



Figure [5] shows the near-infrared spectra of the new 
brown dwarfs. Also plotted for comparison purposes 
is our IRTF/SpeX spectrum of UGPS 0722-05, the 
latest type brown dwarf known previous to this work. 
All of the spectra exhibit deep H 2 and CH4 ab- 
sorption bands characteristic of late-type T dwarfs but 
the J-band peaks of the WISE brown dwarfs are nar- 
rower than the corresponding peak in the spectrum 
of UGPS 0722—05. This peak becomes progressively 
narrower beyond T8 (IWarren et al.ll2007l: iDelorme et all 
I2008al iBurningham et al.ll2008l : iLucas et al.ll2010f >. indi- 
cating that all of the WISE brown dwarfs have spec- 
tral types later than UGPS 0722-05. The spectrum of 
WISEP J1828+2650 is markedly different than that of 
UGPS 0722—05 so we discuss this object in more detail 
in the following section before discussing the other six 
dwarfs. 

4.1.1. WISEP J 1828+2650: The Archetypal Y Dwarf 

The lower panel of Figure [4] shows the 1.15—1.70 
fj,m spectrum of WISEP J 1828+2650 along with the 
spectrum of UGPS 0722-05. The spectrum of 
WISEP J1828+2650, while dominated by the same CH 4 
and H2O absorption bands present in T dwarf spectra, 
has a feature not seen in any T dwarf: the J- and H- 
band peaks, when plotted in units of f\, are essentially 
the same height. As discussed in £ 13.2.41 the spectrum 
of WISEP J1828+2650 is contaminated by light from 
nearby stars. This contamination, which is not removed 
by the aXe software, becomes progressively worse at 
shorter wavelengths (see Figure [3]) such that the true 
spectrum will have an even more extreme J to H band 
peak flux ratio. 

The roughly equal-intensity J and H flux peaks are 
also confirmed by our ground-based, near-infrared pho- 
tometry, which gives J — H = 0.72+0.42 mag. Model 
atmospheres of cool brown dwarfs predict that the near- 
infrared colors, which are blue for the hotter T dwarfs, 
turn back to the red at effective temperatures between 
300 and 400 K as the Wien tail of the spectral energy dis- 
tribution collapses. This turn to the red was proposed 
as one of the tri ggers that might force the creation of a 
Y spectral class (jBurrows et alll2003l !Kirkpatrickll2008lh 

Further underscoring the extreme nature 
of WISEP J1828+2650 is its J-W2 color of 
9.29+0.35 which is over 2 mags redder than the 
WISEP J1541-2250, the second red dest brown dwarf 
in ou r sample at J— W2=7. 18+0.38 (|Kirkpatrick et all 
[20ll . WISEP J1828+2650 is also the reddest brown 
dwarf in our sample in H—W2, J— [4.5], and H— [4.5], 
where [4.5] represe nts the Spitzer Space Telescope 
(IWerner et all 120041) Infrared Array C amera (IRAC; 
iFazio et all 120041 ) channel 2 magnitude. Given the 
extreme nature of both its near-infrared spec- 
trum and near- to mid-infrared colors, we identify 
WISEP J1828+2650 as the archetypal member of the Y 
spectral class. 

4.1.2. WISEP J1738+2732 & WISEPC J1405+ 5534 and 
the 1.5 \xvn NH3 Band 

The 1.15-1.70 (im spectra of WISEPC J1405+5534 
and WISEP J1738+2732 are very similar and yet both 
are distinct from UGPS 0722—05, albeit in less extreme 



The Y dwarfs 



5 



ways than the spectrum of WISEP J1828+2650 (see mid- 
dle and upper panels of Figure [4]). Although the relative 
heights of the J- and ii-band peaks are similar to those 
of UGPS 0722—05, their widths are narrower. The nar- 
rowing of the ii-band peaks is asymmetric, however, as 
most of the change is a result of enhanced absorption on 
the blue wings from 1.51 — 1.58 ^m. What is the under- 
lying cause of this absorption? 

The ii-band spectra of T dwarfs are shaped by CH4 
(and to a lesser extent H2O) longward of 1.6 /xm and 
by H 2 at wavelengths shortward of 1.6 fim. As the 
effective temperature falls, the opacity of the near in- 
frared overtone and combination bands of NH3 becomes 
important since NH3/N2 > 1 for T < 700 K at P=l bar 
(ILodders fe Feglevll20"o"2T) . The emergence of these NH3 
bands has long been suggested as the trigger for a new 
spectral class (IBurrows et all 120031: IKirkpatrickl 120081: 
ILeggett et al.l2007a[) but identifying them has proven dif- 
ficult because they overlap with the strong H2O bands 
and because the abundance of NH 3 can be reduced by 
an order of magnitude due t o vertical mixing in the at- 
mospheres of brown dwa rfs (jSaumon et~alll200l [2006; 
iHubenv fc BurrowslfeOOTD . 

Figure [5] shows the ii-band spectra of the T4, 
T6, and T8 spectral standards, UGPS 0722-05, and 
WISEP J1738+2732 a s well a s the opacities for H 2 Q 
(iFreedman et all [20081) . NH^ dYurchenko et all l20Tlj ), 
and CHa ( |Freedman et all 120081) at T=600 K and P=l 
bar generated by one of us by (R.S.F.). With increas- 
ing spectral type, the blue-wing of the if -band peak be- 
comes progressively suppressed. iDelorme et al.l (|2008a[) 
tentatively identified NH3 absorption on the blue wing of 
the iJ-band spectrum of CFBDS J0059-0114, a T dwarf 
with a spectral type earlier than UGPS 0722—05. How- 
ever, the change in the shape of the blue wing of the H- 
band peak from T6 to UGPS 0722—05 appears smooth, 
suggesting a common absorber or set of absorbers. It 
seems unlikely that NH3 dominates given that it has not 
been identified in the spectra of mid-type T dwarf (T e ff 
~ 1200 K). A sim i lar co nclusion to ours is reached by 
IBurningham et al.l (|2010f ) using spectral indices. 

In contrast, the ii-band spectrum of 
WISEP J1738+2732 stands out in the sequence in 
that it exhibits additional absorption from 1.53 to 1.58 
/im. This absorption broadly matches the position of 
the v\ + v% absorption band of NH3 centered at 1.49 
/jm suggesting that NH3 is the cause of this absorption. 
However, we cannot conclusively identify NH3 as the 
carrier given the low spectral resolution of the data and 
the fact that the absorption lies on the steep wing of the 
H2O band. For example, wate r ice also has an absorp - 
tion band centered at ~1.5 /xm (|Warren fc Brandtll2008l) 
that could potentially produce such absorption if the 
abundance of water ice is high enough. One potential 
avenue for confirming that NH3 is indeed the carrier 
would be to acquire higher spectral reso lution data to 
search for individual NH 3 features (e.g., iSaumon et al.1 
120001: IWarren et al.ll200l . 

4.1.3. The T Dwarf/Y Dwarf Transition 

With WISEP ,11828+2650 classified as the prototyp- 
ical Y dwarf, we can now investigate the transition 
between the T and Y spectral classes. T dwarfs 
are classified at near-infrared wavelengths using the 



iBurgasser et al.l (|2006f) scheme, wherein nine T dwarf 
spectral standards with subtypes rangin g from TO to T8 
are used for direct spectral comparisons. IBurgasser et al.l 
also defined five spectral indices that measure the depths 
of the CH4 and H2O bands which can be used as a proxy 
for direct comparisons. With the discovery of brown 
dwarfs with spectral types later than T8, the question 
of how to extend the IBurgasser et al.l scheme beyond T8 
naturally arises. 

The first >T8 dwarf to be identified was 
ULAS J003402. 77- 005206. 7 (ULAS J0034-0052; 

IWarren et al.l I2007D . Based on a direct comparison 
to the spectrum o f the T8 spect ral standard and 
the values o f the IBurgasser et al.l spectral indices, 
IWarren et al.1 adopted a spectral type of T8.5. A 
second >T8 dwarf soon followed with t he disc overy of 
CFBDS J0059-0114 by IDelorme et all (|2008a[ ). They 
used the Wj index w hich measures the half-width 
of the J band peak (|Warren et alj I2007D and the 
NHs-iJ index which measures the half-width of the 
if-band peak (i.e. the depth of the putative NH3 
absorption), to classify both CFBD S J0059-0114 and 
ULAS J0034-0052 as T9 dwarfs. IBurningham et all 
(|2008f ) added two T dwarfs to the tally of >T8 dwarfs 
with the discovery of ULAS J133553.45+113005.2 
(ULAS J1335+H30) and ULAS J123828.51+095351.3 
(ULAS J1238+0953). Using both direct spectral 
comparison and spectral indice s, they classified t hem 
as T8.5 and T9, respec tively. | Burningham et al.l also 
proposed extending the IBurgasser et al.l scheme to T9 
by assigning ULAS J1335+1130 as the T9 spectral 
standard. Additional T dwarfs with spectral types later 
than T8 have since been discovered (see Table [5]), but to 
date, the latest-type T dwar f currently known is UGPS 
0722—05 (jLucas et al.l 120101 ) which has been classified 
as T10 via a combination of spectral indices and direct 
comparisons to the T9 dwarfs. 

Figure [S] shows the 1.15—1.70 fim spectra of the 
T6, T7, T8 spectral standards, UGPS 0722-05, and 
WISEP J1738+2732. The spectra show smooth changes 
in their spectral morphology with increasing spectral 
type including progressively deeper absorption bands 
centered at 1.15, 1.45, and 1.65 [im and progressively nar- 
rower J- and ii-band peaks. However UGPS 0722—05 
does not appear to be two subtypes later than T8 as 
required by its T10 spectral type. Rather, the changes 
in spectral morphology from T6 to UGPS 0722—05 sug- 
gest that UGPS 0722—05 is more naturally classified as 
a T9. Given its brightness (J=16.5, 1.5 mag brighter 
than ULAS J1335+1130) and position near the celestial 
equator, UGPS 0722—05 also makes an ideal spectral 
standard. We therefore define it to be the T9 spectral 
standard. 

WISEP J1738+2732 is clearly of later type than UGPS 
0722—05 but should it be classified as a T dwarf 
or a Y dwarf? As noted in the previous section, 
WISEP J1738+2732 exhibits excess absorption from 
1.53 to 1.58 /im that we have tentatively ascribed to NH3. 
This absorption becomes even more apparent when the 
spectrum is placed in sequence with the T6 to T9 spec- 
tral standards (lower right panel of Figure 0. Given 
the smooth change in width of the J-band peak and the 
rapid fall in the flux of the blue wing of the ii-band be- 
tween UGPS 0722-05 and WISEP J1738+2732 (which 



6 



Cushing et al. 



suggests the emergence of a new absorption band), we 
classify WISEP J1738+2732 as a Y dwarf and assign it 
a spectral type of YO. Additionally, we tentatively iden- 
tify it as YO spectral standard. 

4.1.4. Classification of the Other WISE Discoveries 

With the T9 and YO spectral standards defined, we 
can return to the question of classifying the other 
new WISE discoveries. The J- and i7-band peaks of 
WISEPC J0148-7202 are slightly narrower than UGPS 
0722-05 and slightly wider than WISEP J1738+2732 
so we classify this dwarf as T9.5. The spectrum 
of WISEPC J1405+5534 is very similar to that of 
WISEP J1738+2732 (see Figure gj). However, we note 
that the wavelength at which the peak iJ-band flux is 
reached is shifted ^60 A to the red relative to UGPS 
0722—05 (see Figure[4]) and the other late-type T dwarfs 
which suggests that WISEPC J1405+5534 may be pecu- 
liar. We therefore classify it as Y0 (pec?). Interestingly, 
a similar, albeit larger, shift of 200 A is seen in the spec- 
trum of Jupiter. 

The spectra of the remaining brown dwarfs, 
WISEPC J2056+1459, WISEP J0410+1502, and 
WISEP J1541-2250, do not have sufficient S/N to con- 
vincingly show the excess absorption from 1.53—1.58 fmx. 
However, the J-band peaks of these three brown dwarfs 
are all narrower than UGPS 0722—05. Indeed the spectra 
of all of them are a better match to the spectral mor- 
phology of WISEP J1738+2732 than UGPS 0722-05 
so we classify these brown dwarfs as Y0 as well. In 
addition, the peak Y-band fluxes of WISEP J0410+1502 
and WISEP J1541-2250 are slightly higher than in 
the spectrum of UGPS 0722—05. This is consistent 
with the bluew ard trend in the Y — J color of late- 
type T dwa rfs (iLeeeett et al.l 120101: IB urningham et aLl 
120101 ) which iBurningham et al.1 ( 2010l ) ascribed to the 
weakening of the strong resonance K I doublet (7665, 
7699 A) as K condenses into KC1. Finally, since the 
spectrum of WISEP J1828+2650 is distinct from both 
UGPS 0722-05 and WISEP J1738+2732, we classify 
it as >Y0. A more precise subtype will require the 
discovery of additional Y dwarfs to bridge the gap in 
spectral morphology between WISEP J1738+2732 and 
WISEP J1828+2650. 

4.1.5. Reclassification of Previously Published >T8.5 Dwarfs 

There are also twelve T dwarfs with spectral types later 
than T8 currently in the literature (see Table Since 
we have reclassified UGPS 0722—05 as a T9 dwarf, wc 
must also reclassify the other eleven dwarfs using this 
new system. To accomplish this, we have smoothed 
the published spectra to a resolving power of i?=150 
and rcsamplcd them onto the same wavelength scale 
as the IRTF/SpeX spectrum of UGPS 0722-05. This 
ensures that differences in resolving power and wave- 
length sampling between the late-type T dwarfs and the 
iBurgasser et al.l IRTF /SpcX spectra of the T dwarf spec- 
tral standards do not adversely affect our classification. 
Table [S] gives the revised spectral types derived from di- 
rect comparison for the twelve T dwarfs with published 
spectral types later than T8. When the J and H band 
regions gave conflicting spectral types, we used the typed 
inferred from the J-band. 



4.1.6. Spectral Indices 

Although the primary (and preferred) method of as- 
signing a spectral type is to compare the spectrum of a 
brown dwarf against that of the spectral standards, the 
use of spectral indices remains popular in the literature. 
We have the refore computed the H 2 Q-J, CH4-J, H 2 Q- 
H, C H4-/J (IBurgasser et al.l [2006]) W T dWarren et al.l 
120071 ). and NH3-/J ijDelorme et alj|2008a|) indices of the 
new WISE brown dwarfs and as well as UGPS 0722-05. 
Figure [7] illustrates the positions of the indices' flux win- 
dows relative to the spectrum of UGPS 0722-05. The 
index values are computed in a Monte Carlo fashion 
whereby 5000 realizations of each spectrum are gener- 
ated by randomly drawing from normal distributions 
with means given by the flux densities at each wavelength 
and standard deviations given by the uncertainty in the 
flux densities. The values of the indices and their uncer- 
tainties are given by the mean and standard deviation of 
the distribution of index values computed from the 5000 
realizations and are listed in Table [4j 

Figure [8] shows the values of the six indices as a func- 
tion of spectral type. Also shown are the index values 
of the T6— T8 spectral standards, a sample of T5— T8 
dwarfs from the SpeX Prism Spectral Library, and the 
twelve T dwarfs with previously published spectral types 
later than T8. The classification of WISEP J1738+2732 
as a Y dwarf is bolstered by the distinctive break in the 
trend of the NH3-i7 values with spectral type which sug- 
gests that a new absorption band has indeed emerged at 
the T/Y dwarf transition. The remaining spectral in- 
dices do not show a break at the T/Y transition, but 
the CH4-J, H2O-H, and Wj indices do show a smooth 
trend with spectral type down to Y0 indicating that they 
can still be used as proxies for direct spectral compar- 
isons. Indeed the value of the Wj spectral index for 
WISEP J1738+2 732 is far from satu r ated s o we support 
the suggestion by IBurningham et all (|2008l) that this in- 
dex can be used as a proxy for direct comparison for 
late-type T dwarfs and early-type Y dwarfs. However, 
the GK4-H index is clearly beginning to saturate at T9 
and the H2O-J index may even reverse at Y0 rendering 
these indices less useful for classification purposes. Fi- 
nally, although there is scatter due to the very low S/N 
of some of the spectra, the new WISE brown dwarfs are 
clearly distinct from the previous >T8.5 dwarfs and clus- 
ter around the Y0 spectral standard. 

4.2. Atmospheric and Structural Properties 
4.2.1. Atmospheric Properties 

In order to investigate the atmospheric properties (e.g., 
T e ff, log g) of the brown dwarfs, we have compared their 
near-infrared spectra to a new preliminary grid of model 
spectra generated with the model atmospheres of Marley 
& Saumon. A detailed descr i ption o f the basic models 
can be found inlMarley et al.l (I2002D . iSaumon fc Marlevl 
(|2^0^ . ICushing et al.1 <\200K and lStephens et al.l (120091) . 
This preliminary grid includes a new NH3 line list 
(jYurchenko et al.l 120111 ) and a new prescription for the 
collision induced opacity for H 2 (Saumon et al., in prep). 
A more detailed study that compares the model spec- 
tra to the near-infrared spectra, and WISE and Spitzer 
photometry is in preparation. 

The grid consists of solar metallicity, cloudless mod- 
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els with the following parameters: =200— 1000 K in 
steps of 50 K; log g = 3.75— 5 in steps of 0.25 (cm 



and K zz =0, 10 cm. s 



Although the opac- 



ities of the condensate clouds are not included in the 
atmospheric models, i.e. they are cloudless, the effects 
on the atmospheric chemistry due to the rainout of the 
condensates is accounted for in the models. This as- 
sumption is reasonable for the silicate and liquid iron 
clouds since they for m well below the obse rvable photo- 
sphere (see however, iBurgasser et al.ll2010l) but may not 
be valid if, as expected, H2O clouds form high in the 
atmosphere of the coldest models. The eddy diffusion 
coefficient, K zz , parametrizes the vigor of mixing in the 
radiative layers of the atmosphere. A value of K zz >0 cm 2 
s _1 results in mixing that can prevent the abundances 
of CO and CH4 (the dominant carbon-bearing species) 
from coming into chemical equilibrium because the mix- 
ing time scales become shorter than the timescales of 
key chem ical reactions involved in the conversion of CO 
to CHa (lLodders fc Ferievl [200l ISaumon et all l20t)l 
iHubenv fc Burrowsll2007l ). Typical values of K zz in the 
stratospheres of gian t planets are 10 2 to 10 5 cm 2 s _1 
(|Saumon et all 12006( 1 The abundances of N 2 and NH 3 
(the dominant nitrogen-bearing species) are also kept 
from coming into chemical equilibrium by mixing, but 
in this case the mixing timescales are set in the convec- 
tive layers of the atmosphere by the mixing length theory. 
As a result, the final non-equilibrium abundances of N2 
and NH3 are not sensitive to variations in the eddy diffu- 
sion coefficient K zz . However by convention, models with 
K zz =0 cm 2 s _1 are in full chemical equilibrium (i.e. the 
effect of convective mixing on the nitrogen chemistry is 
ignored) and models with K zz ^ cm 2 s _1 exhibit both 
carbon and nitrogen non-equilibrium chemistry. 
The best fitting models are identified using the 

goodness-of-fit statistic, Gk = J27=i Wi ■ 
where n is the number of data pixels; Wi is the weight for 
the ith wavelength (set to unity in this case); ft and Fkj 
are the flux densities of the data and model fc, respec- 
tively; <7j; are the errors in the observed flux densities; 
and Cfc is an unknown multiplicative constant equal to 
(R/d) 2 , where R is the radius of the st ar and d is the 
distance to the star ([Cushing et a l. 2008) . In order to in- 
crease the S/N of the data, we first smoothed the higher 
resolution spectra to i?=200. The model spectra were 
also smoothed to the same resolving power and linearly 
interpolated onto the wavelength scale of the data. For 
each model we compute the scale factor Ck by minimiz- 
ing Gk with respect to Ck and identify the best fitting 
model as having the global minimum Gk value. To esti- 
mate the range of models that fits the data well, we run 
a Monte Carlo simulation that uses the uncertainties in 
the individual spectral points and the uncertainties in the 
absolute flux calibration of the spectra to generate 10 4 
simulated noisy spectra. The fitting process is repeated 
on each simulated spectrum and models that are consis- 
tent with the best fitting model at the 3<r level are con- 
sidered equally good representations of the data. We did 
not attempt to fit the spectrum of WISEP J1828+2650 
because it is contaminated with light from other stars in 
the WFC3 field of view (see 33T4l and EXT]). After dis- 
cussing the results of the fits to the spectra of the other 



brown dwarfs, we return to estimate an approximate ef- 
fective temperature for WISEP J1828+2650. 

Table [6] lists the best fitting model parameters for each 
brown dwarf, as well as UGPS 0722—05. The derived 
effective temperatures of the WISE brown dwarfs are all 
cold, ranging from 350 to 500 K. Indeed all but one have 
estimated effective temperatures of less than 450 K mak- 
ing them the coldest spectroscopically confirmed brown 
dwarfs known. Five out of the six best fitting models 
also have K zz ^0 which indicates that vertical mixing is 
present in the atmospheres of these cold brown dwarfs. 
This is not a surprising result given that strong evidence 
for vertical mixing in the atmospheres of brown dwarfs 
has been found at longer wave l engths (ISaumon et all 
2001 [20071: iLeggett et al]|2007bl IBurgasser et all 120081: 



Leggett et all 1200ft iStephens et all 120091: iGeballe et al l 
200ft iCushing et all l20l6T T At such low temperatures, 



the effects of non-equilibrium chemistry on the J- and 
-ff-band spectra of brown dwarfs is limited to weaken- 
ing the NH3 absorption bands. The detection of mixing 
therefore underscores the fact that NH3 is probably at 
least partially responsible for the absorption seen on the 
blue wing of the ff-band peak of WISEPC J1405+5534 
and WISEP J1738+2732. 

Figure [H] shows the best fitting model spectra overplot- 
ted on the data. Since this is the first time such cold 
model spectra have been compared to observed spectra, 
the agreement between the models and the data is en- 
couraging. In particular, the height and width of the 
J-band peaks are well matched by the model spectra. 
Previous studies of hotter brown dwarfs fail to match 
both the peak and width of this peak (e.g. . ILeggett et ail 
1200ft IBurgasser et al.ll20lil ). The improved fits may be a 
result of the fact that we are fitting a smaller wa velength 
range (jCushing et all 120081: fSeifahrt et all 1200ft ) and/or 
because the high J vibration-rotation lines (the so-called 
"hot" lines), whose cross-sections are less well known, 
become less important at such cold temperatures. 

The models do, however, provide a poor fit to the 
blue wing of the .ff-band peak of the spectra and 
fail to reproduce the heights of the F-band peaks of 
WISEP J0410+1502 and WISEP J1541-2250. Note 
that the peak of the Y-band is shaped by the 2vi+2v<± 
band of NH3 centered at about 1.03 /zm and therefore Y- 
band spectra of cold brown dwarfs could provide the first 
clear detection of NH3 at near-infrared wavelengths. In 
principle, the blue wing of the .ff-band model spectrum 
could be brought into better agreement with the data 
by further reducing the abundance of NH3. However, 
as noted above, the abundance of NH3 is insensitive to 
variations in K zz because it is quenched in the convective 
region where the mixing time scale is set by the mixing 
length theory and not by the eddy diffusion coefficient. 
Therefore, the mismatch between the data and models 
is most likely a result of some other inadequacy in the 
model atmospheres. 

Finally, although we cannot fit the models to the spec- 
trum of WISEP J1828+2650, we can still estimate a 
rough effective temperature. The most salient feature 
of the spectrum is that the J- and .ff-band peaks are 
roughly the same height in flux density units of f\. 
Only model spectra with T c & < 250 K have J-band 
peak fluxes that are equal to or less than the ff-band 
peak fluxes. A second estimate of the effective tern- 



Cushing et al. 



perature can be derived from the observed J— W2 color 
of 9.29±0.35. We computed synthetic Mauna Kea Ob- 
serva tories Near-Infrared (MKO-NIR, iTokunaga et al.l 
120021 ) J and W2 magnitudes for each model in the 
grid and find that model spectra with T e fj=275— 300 
K have J— W2 colors that fall within ±2er of the ob- 
served color. Taken together, these estimates suggest 
that an appropriate upper limit to the effective temper- 
ature of WISEP J1828+2650 is ~300 K which makes 
WISEP J1828+2650 the coolest spectroscopically con- 
firmed brown dwarf known. 

4.2.2. Structural Properties 

With estimates of the effective temperatures and sur- 
face gravities of the new brown dwarfs in hand, we can 
also estimate their radii (R) and masses (M) using evo- 
lutio nary models. We used the cloudless structure mod- 
els of lSaumon fc Marlevi (|2008ft because they used atmo- 
spheric models that are nearly identical to the ones we 
used in our analysis for boundary conditions. As a re- 
sult, the derived T e g, log g, R, and M estimates are all 
self-consistent. The radii and masses of the brown dwarfs 
are given in Table 

4.3. Spectroscopic Distance Estimates 

The value of the multiplicative constant Ck~(R/d) 2 
derived as a by product of the atmospheric model fit- 
ting procedure can be used to estimate a so-called 
"spectroscopic distance" (rf S nsr.) to brown dwarfs if their 
radii can be determined (e.g., iBowler et al.l [2009h . In 
the absence of direct measurements of brown dwarf 
radii, we can use the radii computed using evolution- 
ary models and (T e g, log g) values. The spectroscopic 
distances of the new WISE brown dwarfs and UGPS 
0722-05 are given in Table Also listed in Table 
are the photome tric distance estimates o f the WISE 
brown dwarfs from iKirkpatrick et all (120111) and paral- 
lactic distances to UGPS 0722-05 dLucas et al. l 12010ft 
and WISEP J1541-2250 (IKirkpatrick et all 120 lift . The 
former distance estimates are computed using a spectral 
type-W2 relation derived from known brown dwarfs with 
spectral types ranging from L0 to T9 and with 7r/ov > 
3. The photometric distances of the new WISE brown 
dwarfs are based on an extrapolation of this relation and 
therefore should be treated with caution. 

The agreement between the three distance estimates 
range from good to poor. For example, the spectro- 
scopic and parallactic distances of WISEP J 154 1—2250 
are in good agreement but the photometric distance is 
discrepant by a factor of two to four. Perhaps most 
discouraging is the mismatch between the spe ctroscopic 
and p arallactic distances of UGPS 0722—05. iLiu et all 
([2011ft recently showed that the agreement between the 
spectroscopic distances (derived using only near-infrared 
spectra) and the parallactic distances of ten late-type T 
dwarfs range from 10% to a factor of two, with no appar- 
ent trend with spectral type or distance. This suggests 
that spectroscopic distances should only be used to con- 
firm that the physical properties of brown dwarfs (T e g, 
R, M) derived from atmospheric and evolutionary mod- 
els are consistent with the known distance to the brown 
dwarf. 

A corollary to this statement is that if the spectroscopic 
and parallactic distances arc discrepant then some com- 



bination of the T e ff, logg and R values are in error. In 
order to estimate the significance of the bias in the spec- 
troscopic distance estimate introduced by systematic er- 
rors in the inferred atmospheric properties, we have run 
a 1—2.5 (im model with T e ff=600 K, log</=4.5 cm s~ 2 
through the fitting procedure described in §4.2.11 The 
model spectrum was first multiplied by an appropriate 
value of {R/d) 2 corresponding to 10 pc. Figure [TOl shows 
the ratio of d S pcc/ 10 pc derived for models with effec- 
tive temperatures from 500 to 700 K and surface grav- 
ities from 3.75 to 5.0. The maximum change in d spcc 
for these models is approximately a factor of two for a 
change of +100 K/-0.75 dex and -100 K/+0.5 dex in 
Teff/logg. The spectroscopic dista nce is also most sensi- 
tive to changes in T e g as noted by ILiu et all ([201 lft . 

The apparent mismatch between the spectroscopic and 
parall actic distance estima tes is perhaps not so surpris- 
ing as lCushing et al.l ([2008ft have shown that variations of 
order 100 to 200 K are common when estimating the ef- 
fective temperatures of L and early- to mid- type T dwarfs 
using spectra that cover only a fraction of the spectral 
energy distribution. These variations are most likely ex- 
acerbated by the fact that only ^35% (in flux density 
units of f^) of the bolometric flux of a T o fj=600 K brown 
dwarf is emitted at near-infrared wavelengths. In sum- 
mary, given the uncertainties in the model atmospheres 
and the difficulty in estimating the effective temperatures 
and surface gravities of brown dwarfs, it is not surpris- 
ing that spectroscopic distance estimates do not always 
agree with the parallactic distances. 

5. DISCUSSION 

The new WISE brown dwarfs presented herein are the 
coldest (T G ff=300— 500 K) spectroscopically confirmed 
brown dwarfs currently known. However as noted in §1, 
WD 0806-661B and CFBDSIR J1458+1013B have esti- 
mated effective te mperatures of ~300— 400 K based on 
photometry alone dLuhman et all 120111 : ILiu et ail 120111 : 
iRodriguez et all 120111 ). How do the properties of these 
two brown dwarfs compare with the new WISE brown 
dwarfs? 

The upper panel of Figure [Tl] shows the absolute J- 
band magnitude (Mj) as a func tion of spectral typ e 
for a sample of field T dwarfs (|Leggett et all 12010ft . 
WISEP J1541-2250 (the only WISE brown dwarf with 
a measured parallax), WD 0806— 661B (L uhman et 
al. s ubmitted) and CFBDSIR J1458+1013B ([Liu et all 
12011 . The value of Mj increases precipitously be- 
yond a spectral type of T8 and reaches ~23.9 for 
WISEP J1541-2250 (Y0). The absolute magnitude of 
CFBDSIR J1458+1013B falls between the two T9 dwarfs 
and WISEP J1541— 2250 suggesting that it has a spec- 
tral type of T9— Y0. However the absolute magnitudes of 
more late-type T dwarfs and Y dwarfs must be measured 
before any strong conclusions can be drawn based on ab- 
solute magnitudes alone. The absolute magnitude of WD 
0806-661B is still only a limit (Mj > 22.5) which leaves 
open the possibility that WD 0806— 661B is even fainter, 
and thus presumably cooler than, WISEP J1541— 2250. 
Either way, it is clear that based on the trend of Mj 
with spectral type that observing even colder objects at 
near- infrared wavelengths is going to become increasingly 
difficult unless they are very close the Sun. 

The lower panel of Figure [Tl] shows the J — H col- 
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ors as a fu nction of spectral t ype for a sample of field 
T dwarfs (|Leggett et all l2010f ). the new WISE brown 
dwarfs, and CFBDSIR J1458+1013B (WD 0806-661B 
has yet to be detected in either the J or H band). 
Some of the colors of the WISE brown dwarfs have 
large uncertainties so we also comp uted synthetic MKO - 
NIR J — H colors as described in iRavner et al.l (|2009f ); 
they are shown as triangles in Figure [11] Since the 
WFC3 / HST spectra do not span the entire H band, we 
used the spectrum of UGPS 0722—05 to extend the spec- 
tra of WISEP J1738+2732 and WISEPC J1405+5534to 
the limit of the H band filter. The synthetic and ob- 
served photometry of WISEPC J1405+5534 are clearly 
discrepant and it is unclear what the underlying cause 
is. The scatter in the J — H colors of the dwarfs at 
the T/Y transition makes it difficult to assign a spectral 
type to CFBDSIR J1458+1013B but it is broadly consis- 
tent with T6— Y0. The overall trend with spectral type 
suggests suggests that the J — H colors may be turn- 
ing towards to the red at the T/Y transition. This turn 
is broadly consistent with theoretical models which pre- 
dict that the J — K color also turns to wards the red at 
r efi =390-450 K ([Burrows et al.ll200l . However given 
the small number of objects and the large uncertainties 
in colors, a definitive conclusion can not yet be reached. 

Finally, given the rapid increase in the absolute J- 
band magnitude at the T/Y transition, it is reason- 
able to ask whether the near-infrared is the appropri- 
ate wavelength range with which to define the Y spec- 
tral class. Indeed historically as cooler and cooler stars 
were discovered, the wavelength range used for spectral 
clas sification moved from the blue violet at 3900—4900 
A ([Morgan et al.l 119431) . through the red optical at 
7000-10000 A (jBoeshaarjll9"7l iKirkpatrick et~aT1 Il99l"[ 
11999ft. and finally into the near infrared at 1—2.5 /an 
([Burgasser et al.l 120061) . Since Y dwarfs emit the ma- 
jority of their flux in the mid-infrared, it seems only 
natural to devise the spectral classification system at 
these wavelengths. The smooth spectral morphological 
changes seen i n the 5.5—14.5 /xm Spitzer spectra of L 
and T dwarfs ([Cushing et al.ll2006h suggest that a mid- 
infrared classification scheme for the Y dwarfs is plau- 
sible. Unfortunately, observing at wavelengths longer 
about 2.5 /im from the ground is exceedingly difficult 
due to the high thermal background. Observations con- 
ducted from space do not suffer from this limitation but 
there are currently no space-based facilities (including 
Spitzer) capable of obtaining mid- infrared spectra of cold 
brown dwarfs. The James Webb Space Telescope will 
provide a platform with which to observe cold br own 
dwarfs (|Burrows et al.ll2003t Parley fc Leggetti2009[) but 
its launch is, at best, still years away. 

We are therefore left in the unfortunate position of ei- 
ther waiting for a (space- or ground-based) facility capa- 
ble of sensitive mid-infrared observations or constructing 
a classification scheme in the near infrared. Given the 
large number of cold brown dwarfs now known, we be- 
lieve it is important to devise a system with which to 
classify them based on the data currently available. The 
creation of a near-infrared scheme in no way invalidates 
any future mid-infrared system that may be devised. In- 
deed the classification of brown dwarfs at two different 
wavelengths is not unprecedented as both the L and the 



T dw arfs have classification systems based in the re d op- 
tical ([Kirkpatrick et al.lll999HBurgasser et al.ll2003ft and 



the near infrared (iGeballe et al.l 20021 : IKirkpatrick et al.l 
120101: IBurgasser et all 120061 ). Ultimately, the utility of 
any classification system will be measured by whether or 
not it is adopted by the brown dwarf community. 

6. SUMMARY 

As part of an ongoing search for the coldest brown 
dwarfs in the solar neighborhood using the Wide-field 
Infrared Survey Explorer, we have discovered seven ul- 
tracool brown dwarfs whose near-infrared spectra indi- 
cate that they are the latest-type brown dwarfs currently 
known. Based on the spectral morphological differences 
between these brown dwarfs and the late- type T dwarfs, 
we have identified six of them as the first members of 
the Y spectral class. A comparison to the model spec- 
tra of Marley & Saumon indicate that they have effec- 
tive temperatures ranging from 300 to 500 K and thus 
are the coolest spcctroscopically confirmed brown dwarfs 
currently known. 
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FIG. 1.— 2x2 arcmin DSS I, 2MASS J and H, WISE Wl, W2, and W3, and a W1W2W3 false color composite of the five new WISE 
brown dwarfs. In the color composite images on the far right, the Wl, W2, and W3 are bands are color coded blue, green, and red, 
respectively. 
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FlG. 2. — Near-infrared spectra of the new WISE brown dwarfs (black) as compared to the spectrum of UGPS 0722—05 (red). The data 
have been normalized to unity at the peak of the J band (except for WISEP J1828+2650 which is normalized to unity at the peak of the 
H-band) and offset by constants (dotted lines). 
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FIG. 3.— Top: Subimage of the drizzled WFC3/HST grism image centered on the position of WISEP J1828+2650. The location of the 
spectrum of WISEP J1828+2650 is indicated in red along with the positions of the J- and H-band peaks. The location of the contaminating 
starlight is shown in purple and consists of second and third order light from two other stars in the WFC3 field of view. Bottom: The 
spectrum of WISEP J1828+2650 (red) and the contamination spectrum (purple). The stellar contamination becomes progressively worse 
at shorter wavelengths. 
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FIG. 4.— 1.15-1.70 /im spectra of WISEP J1738+2732, WISEPC J1405+5534, and WISEP J1828+2650 along with the spectrum UGPS 
0722—05. The uncertainties in the spectra are shown as gray bars. The spectra were all normalized to unity at the peak of the H-band 
(1.58 (im). Prominent molecular absorption bands are indicated. 
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Fig. 5.— Top: _ff-band spectr um of 2MASS J2 2541892+3123498, SDSS J162414.37+002915.6 and 2MASS J04151954-0935066, the T4, 
T6, and T8 spectral standards l IBurgasser et al.ll200rT) . along with the spectrum of UGPS 0722-05, and WISEP J1738+2 732. The spectra 
have been no rmalized to unity at th eir peak flux value s. Bottom: Opac ity data computed in chemical equilibrium for NH3 ( Yurchcnko ct al. 
[20TOf) . H 2 HFreedman et al.|[2008T) , and CH 4 llf'Veedman et al.|[2fJ08T) at T=600 K and P=l bar. Note that the change in the spectral 
morphology of the blue wing of the //-band peak is similar between T6/T8 and T8/T9 suggesting a common absorber or set of absorbers. 
In contrast, the spectrum of WISEP J1738+2732 exhibits excess absorption that matches the position of the NH3 absorption shown in the 
lower panel. 
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FIG. 6.— IRTF/SpeX spectra of the IBureasser eFaLl (2006) spectral standards, SDSS J162414.37+002915.6 (T6), 
2MASS J07271824+1710012 (T7), and 2MASS J04151954-0935066 (T8), our IRTF/SpeX spectrum of UGPS 0722-05 and the 
WFC3/ HST spectrum of WISEP J1738+2732. The spectra have been normalized to unity at their peak flux level in each panel. 
Prominent molecular absorption bands are indicated in the top panel. 
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FIG. 7.— Illustration of the H2O-J, Wj, CH4-J, CH4-J, NH3-H, and CH4-H indices overplotted on the spectrum UGPS 0722-05. 
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FIG. 8.— Values of the H 2 0- J, CH 4 - J, H 2 0-H, CH 4 -if IIBurgasser et al.12009) . Wj llWarren et alJ2007l) . and NH 3 -_ff IIDelorme et al.l 
2008a) spectral indices as a function of spectral type. The black points are for the T6— Y0 spectral standards. The grey points were computed 
using spectra of late-type (T5— T8) T dwarfs from the SpeX Prism Spectral Library. The red points are the twelve T dwarfs with previously 
published spectral types later than T8 and the green points are the remaining six WISE brown dwarfs. For plotting purposes only, we have 
assigned WISEP J1828+2650 a spectral type of Yl. 
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FlG. 9. — Best fitting models (red) overplotted on the near-infrared spectra of six of the seven new WISE brown dwarfs (black). The 
spectra were normalized to unity at the peak flux in the J-band and offset by constants (dotted lines). The uncertainties in the spectra 
are given by grey bars. The best fitting model parameters are given in the form T e ff (K)/log g (cm s~ 2 )/K z2 , (cm 2 s" 1 ). 
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TABLE 1 
WISE Photometry 



Object 


Wl 


W2 


W3 


W4 


W1-W2 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


WISEPC J014807.25-720258.8 


18.812±0.529 a 


14.584+0.052 


>12.579 


>9.521 


4.228+0.532 


WISEP J041022.71+150248.5 


>18.101 


14.190+0.059 


12.472+0. 482 a 


>8.923 


>3. 911+0.059 


WISEPC J140518. 40+553421. 5 


>17.989 


14.085+0.041 


12.312+0.252 


>9.115 


>3. 904+0.041 


WISEP J154151. 65-225025. 2 


>17.018 


13.982+0.112 


12. 134+0. 443 a 


>9.064 


>3. 036+0.112 


WISEP J173835. 52+273258.9 


18.155+0.362 


14.535+0.057 


12.536+0.350 


>9.182 


3.620+0.366 


WISEP J182831. 08+265037.8 


>18.452 


14.276+0.050 


12.320+0.291 


9.147+0. 438 a 


>4. 176+0.050 


WISEPC J205628. 90+145953.3 


> 17.742 


13.852+0.043 


11.791+0.222 


>8.646 


>3. 890+0.043 



NOTE. — Objects designated as WISEP arc from the Preliminary Release Source Catalog while objects designated as 
WISEPC arc from the first-pass processing operations coadd Source Working Database. Magnitudes are in the Vega system 
and arc based on profile fits (wlmpro, w2mpro, w3mpro, w4mpro). Upper limits are at the 95% confidence level (see 
http : //wise2 . ipac . caltech. edu/docs/release/prelim/expsup/sec4_3c .html#ul2}. 

a S/N<3 



TABLE 2 
Near-Infrared Photometry 



Object 


Filter 


Instrument 


Magnitude 


Exposure 


Coadds 


Number 


Total 


Date 








(mag) 


Time (sec) 




of Images 


Exp. (sec) 


(UT) 


WISEPC J0148-7202 


J 


PANIC 


18.96+0.07 


30 


1 


18 


540 


2010 Aug 01 




H 


PANIC 


19.22+0.04 


15 


1 


108 


1620 


2010 Aug 01 


WISEP J0410+1502 


J 


WIRC 


19.25+0.05 


60 


1 


15 


900 


2010 Aug 29 




H 


WIRC 


19.05+0.09 


30 


1 


15 


1800 


2010 Jul 26 


WISEPC J1405+5534 


.1 


WIRC 


20.20+0.13 


30 


2 


15 


1800 


2010 Jul 26 




H 


WIRC 


21.45+0.41 








5400 


multiple 


WISEP J1541-2250 


J 


NEWFIRM 


21.16+0.36 


30 


2 


10 


600 


2011 Apr 17 




H 


NEWFIRM 


20.99+0.52 


5 


12 


10 


600 


2011 Apr 17 


WISEP J1738+2732 


J 


WIRC 


19.47+0.08 


60 


1 


15 


900 


2010 Jul 26 




H 


WIRC 


20.66+0.38 


30 


2 


15 


900 


2010 Jul 26 


WISEP J1828+2650 


J 


NIRC2 


23.57+0.35 


120 


1 


6 


720 


2010 Jul 01 




H 


NIRC2 


22.85+0.24 


120 


1 


9 


1080 


2010 Jul 01 


WISEPC J2056+1459 


.1 


WIRC 


19.31+0.12 


60 


1 


15 


900 


2010 Aug 29 




H 


WIRC 


>19.5 


30 


2 


15 


900 


2010 Aug 29 




J 


NIRC2 


19.21+0.07 


120 


1 


9 


1080 


2010 Jul 01 




H 


NIRC2 


19.56+0.18 


120 


1 


6 


720 


2010 Jul 01 



TABLE 3 
Spectroscopy Log 



Object 


Instrument 


UT Date 


Mode 


Slit Width 


Int. Time 


AO V Calibrator 










(arcsec) 


(sec) 


Star 


WISEPC J0148-7202 


FIRE/Magellan 


2010 Sep 18 


Longslit 


0.6 


960 


HD 1881 


WISEP J0410+1502 


FIRE/Magellan 


2010 Nov 18 


Longslit 


1.0 


600 


HD 18620 


UGPS 0722-05 


SpeX/IRTF 


2011 Jan 26 


LowResl5 


0.5 


1440 


HD 50931 


WISEPC J1405+5534 


WFC3/HST 


2011 Mar 14 


G141 




2212 




WISEP J1541-2250 


FIRE/Magellan 


2011 Mar 27 


Longslit 


0.6 


1522 


HD 130755 


WISEP J1738+2732 


WFC3/ffi>T 


2011 May 12 


G141 




2012 




WISEP J1828+2650 


WFC3/ffi>T 


2011 May 09 


G141 




2012 




WISEPC J2056+1459 


NIRSPEC/Kcck 


2010 Oct 21 


low-res (N3) 


0.38 


2400 


HD 198070 




NIRSPEC/Keck 


2010 Nov 22 


low-res (N5) 


0.38 


1800 


HD 198069 



2G 



Cushing et al. 



TABLE 4 
Spectral Indices 



Object 


Spectral 


H2O-J 


CH4-J 


h 2 o-h 


CiU-H 


Wj 


NH 3 -iJ 




Type 















UGPS 0722-05 a T9 +0.009 (0.004) +0.115 (0.003) +0.115 (0.007) +0.075 (0.005) +0.215 (0.003) +0.527 (0.008) 

WISEPC J0148-7202 T9.5 -0.017 (0.011) +0.076 (0.006) +0.070 (0.013) +0.016 (0.012) +0.152 (0.006) +0.431 (0.015) 

WISEP J0410+1502 Y0 -0.043 (0.090) +0.036 (0.059) -0.036 (0.458) -0.105 (0.434) +0.155 (0.052) +0.380 (0.574) 

WISEPC J1405+5534 Y0 (pec?) +0.020 (0.026) +0.031 (0.017) +0.063 (0.025) +0.162 (0.028) +0.099 (0.017) +0.346 (0.027) 

WISEP J1541-2250 Y0 -0.046 (0.115) +0.040 (0.060) +0.155 (0.216) +0.057 (0.208) +0.081 (0.066) +0.259 (0.200) 

WISEP J1738+2732 YO +0.036 (0.008) +0.050 (0.005) +0.045 (0.008) +0.050 (0.009) +0.149 (0.005) +0.349 (0.010) 

WISEP J1828+2650 >Y0 +0.248 (0.152) +0.257 (0.099) +0.049 (0.091) +0.129 (0.101) +0.515 (0.106) +0.503 (0.110) 

WISEPC J2056+1459 YO ■-- +0.030 (0.009) ■■■ +0.050 (0.034) +0.129 (0.012) +0.240 (0.031) 



Note. — The H 2 — J and H 2 0-H 



ot be 



nputcd for WISEPC J205628. 90 + 145953.3 be 



isc its spectr" 



2S not span the entire - 



clcngth : 



nge of the 



The values differ from that measured by Lucas et al. 12010). Our two spectra agree well except deep in the CH 4 and H 2 absorption bands, where our spectrum exhibits 
jwer flux levels. The reason for this discrepancy is unclear but it may be a result of the fact that the Lucas et al. spectrum was created by merging separate spectra that 
•ere absolutely flux calibrated using near-infrared photometry. 



TABLE 5 

Previously Published Brown Dwarfs with Spectral Types Later than T8 



Object 


Previous 


Reference 


Adopted 




Spectral Type 




Spectral Type 



Ross 458C 

ULAS J123828.51+095351.3 
ULAS J130217.21+130851.2 
ULAS J003402. 77-005206.7 

CFBDS J005910.90-011401.3 

WISEPC J045853. 90+643451.9 
UGPS J052127.27+ 364048. 6 
ULAS J133553.45+113005.2 
Wolf 940B 

WISEPC J181210. 85+272144.3 
UGPS J072227.51-054031.2 
CFBDSIR J145829+101343AB 



T8 
T8.5p 
T8.5 
T8.5 
T8.5 

T9 
>T9 

T9 

T9 
T8.5 

T9 
T8.5 
T8.5: 
T10 
T9.5 



Burgasser^e^al^ (2010) 
BurninghanL.e^aL (2011a) 
Burninghani^itjiL (2008) 
13jmriiiglmnirtjiL (2010) 
Warren et al. (2007) 
Burningham et al. (2008) 
Delorme et al. (2008) 
Burninehani^taL (2008) 
Mainzer et al. (2011) 
BmminghametaL (2011b) 
Burningham et al. (2008) 
Burningham et al. (2009) 
BurgasseretaL (2011) 
Lucas et al. (201oT ~ 
Liu et al. (2011) 



T8 



T8.5 

T8.5 

T8.5 

T8.5 

T8.5 

T8.5: 

T9 

T9 



TABLE 6 

Atmospheric and Structural Properties 



Object 


SpType 


T e g 




log 9 




log X zz 


R 




M 






(K) 




(cm s — 


2 ) 


(cm 2 s -1 ) 


(Rjup 


) 


(M Jup ) 


UGPS 0722-05 


T9 


650 




4.00 (4.00- 


-4.25) 


4 


1.21 (1.14- 


-1.21) 


6 (6-9) 


WISEPC J0148-7202 


T9.5 


500 (500- 


-500) 


4.50 (4.50- 


-4.75) 


1 


1.04 (0.96- 


-1.04) 


13 (13-21) 


WISEP J0410+1502 


Y0 


450 (400- 


-500) 


3.75 (3.75- 


-4.25) 





1.22 (1.09- 


-1.22) 


3 (3-9) 


WISEPC J1405+5534 


Y0 (pec?) 


350 




5.00 




1 


0.86 




30 


WISEP J1541-2250 


Y0 


350 




4.50 (4.25- 


-4.5) 


4 


1.01 (1.01- 


-1.07) 


12 (8-12) 


WISEP J1738+2732 


Y0 


350 (350- 


-400) 


4.75 (4.75- 


-5.00) 


4 


0.93 (0.86- 


-0.94) 


20 (20-30) 


WISEP J1828+2650 


>Y0 


<300 
















WISEPC J2056+1459 


Y0 


350 (350- 


-400) 


4.75 (4.50- 


-5.00) 


4 


0.93 (0.86- 


-1.01) 


20 (12-30) 



NOTE. — The parameters for the best fitting Marley 8z Saumon models arc listed and the range of parameters consistent 
with the data is given in parentheses. The effective temperature limit for WISEP J1S28+2650 was estimated by identifying 
those models with peak J-band fluxes equal to or less than the peak flux in the H band and by comparing the observed J-W2 
color to model J— W2 colors. 



The Y dwarfs 



TABLE 7 
Distance Estimates 



Object SpType d spec (pc) a (pc) b d phot (pc) c 



UGPS 0722-05 T9 11.1(10.4-11.1) 3.6-4.7 

WISEPC J0148-7202 T9.5 14.7(13.1-14.7) ■•■ 12.1 

WISEP J0410+1502 Y0 7.1 (3.3-8.7) ■ ■ • 9.0 

WISEPC J1405+5534 Y0 (pec?) 3.8 • • • 8.6 

WISEP J1541-2250 Y0 1.8 (1.8-2.0) 2.2-4.1 8.2 

WISEP J1738+2732 YO 3.4(3.4-7.3) ■•■ 10.5 

WISEP J1828+2650 >Y0 ■ ■ • ■ • ■ <9.4 

WISEPC J2056+1459 YO 3.0 (2.4-6.4) ■ ■ • 7.7 



a Spectroscopic distance estimates derived as described in EI4.2.1I The distance corre- 
sponding to the best fitting model is given and the range of distances corresponding 
to models that are consistent with the data are given in parentheses. 
b Par allact ic distance for UGPS 07 22-05 and WISEP J1541-2250 from lLucas et al.l 

120101) and lKirkpatrick et all 120111) . r espectively. 

c Photometric distance estimates from.Kirkpatrick et al. 1 120111) . 



